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We present a calculation of a neutrino decay scenario 
in the early Universe. The specific decay is 1/2 —* vi + 0, 
where (f} is a, boson. If there is a neutrino mass hierarchy, 
m,^^ < TUv^ < TUv^, we show that it is possible to generate 
stimulated decay and effects similar to atomic lasing without 
invoking new neutrinos, even starting from identical neutrino 
distributions. Under the right circumstances the decay can 
be to very low momentum boson states thereby producing 
something similar to a Bose condensate, with possible conse- 
quences for structure formation. Finally, we argue that this 
type of decay may also be important other places in early 
Universe physics. 

95.35, 13.35.H, 98.62.L 



I. INTRODUCTION 

Any neutrino that is absolutely stable must have a 
mass less than roughly 30 eV in order not to overdose 
the Universe |l|. On the other hand, it is very plausible 
that neutrinos have a non-zero mass. This makes neu- 
trino decays necessary if their mass is larger than the 
cosmological limit for stable neutrinos. There are several 
possible decay modes |y], for example the simple radia- 
tive decay V2 —^ i/i -I- 7 or the flavor changing neutral 
current weak decay V2 -^ Fit^iJ/i. Another possibility is 
the decay 



H ^ F 



(1) 



where H \s a. heavy neutrino, i^ is a light neutrino and 
(^ is a light scalar particle (for example the majoron). It 
is not essential that the boson should have spin 0, but 
it simplifies calculations significantly. All of these decay 
modes have been investigated in an early Universe con- 
text g . We choose to focus on the last decay type in or- 
der to investigate some very interesting non-equilibrium 
effects that pertain to this decay. It turns out that under 
the right circumstances the heavy neutrino will decay by 
stimulated decay into F and (j) instead of just decaying 
thermally. This runaway process leads to the formation 
of a very low momentum component of the scalar par- 
ticle, (j), while the rest are in an almost thermal distri- 
bution. The scenario has been described by Madsen 0] 
and subsequently in more detail by Kaiser, Malaney and 
Starkman [p|,pl, the idea being that the 4> particle could 
make up the dark matter in our Universe. One of the 
currently favoured models of structure formation is the 
so-called Mixed Dark Matter (MDM) model 0. This is 



essentially a Cold Dark Matter (CDM) model, but with 
a small matter content of Hot Dark Matter (HDM) that 
lessen the fluctuations on small scales. The (j) particles 
could fit both roles, the cold component acts as cold dark 
matter and the thermal component as hot dark matter. 

There are, however, several problems with this sce- 
nario. First of all, the favored mixed dark matter mod- 
els consist of roughly 5% baryons, 25% HDM and 70% 
CDM. If there is too much HDM, structure formation 
will be affected so as to produce too little clustering on 
small scales (if normalised to COBE data). The sim- 
ple scenario of Refs. |#-Hl predicts only about 35% CDM 
and 60% HDM, which is incompatible with the observed 
structure [|j. 

We develop further this model and in Section 2 we list 
the relevant equations for calculating the particle distri- 
bution functions resulting from the decay. In Section 3 we 
describe a new way of producing neutrino lasing that is 
very similar to an atomic 3-level laser. Here it is possible 
to achieve lasing without invoking new neutrinos and to 
increase the CDM fraction, which however still remains 
well below 70%. Section 4 contains a description of the 
numerical results from our investigation of this 3-lcvel 
model and, finally. Section 5 contains a discussion of the 
possible consequences of this type of decays as well as a 
discussion of various astrophysical limits on them. From 
this discussion it will become clear that lasing decays are 
most likely to take place after the epoch of nucleosynthe- 
sis (r ~ 1 — O.lMeV). Of course they are also required 
to take place before matter-radiation equality (T ~ leV) 
if the bosons are to act as ordinary MDM. 



II. THE BASIC SCENARIO 

The fundamental equation that describes the evolution 
of all particle species during the decay is the Boltzmann 
equation. 



L[/]-Cdec[/]+Ccou[/]. 



(2) 



The left hand side is the basic Liouville operator in an 
expanding Universe 



L[/] 



9/ 
dt 



dR I df 
HR^dp' 



(3) 



The right hand side of Eq. (E) involves the specific physics 
of the decay process as well as possible other collision 
processes such as elastic scattering and annihilation. We 
shall assume that the term Ccoii can be neglected, but 



possible limitations to this assumption will be discussed 
later. For simplicity we assume that the decay is isotropic 
in the rest-frame (see, however, Refs. ||,g] for a discussion 
of the case of anisotropic decays) . 

The decay terms for the different particles have been 
calculated in for example Refs. |^,||. They are of the 
form 
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dE^A{fHjF,U), (4) 



dEHHiHjF.U), (5) 



dEHA{fHjFj4>), (6) 



where AifnjFj^,) = /ff(l - /f)(1 + U) - IfMI 
Jh), ml^mjj- 2{ml + to|) + (m^ - mD^/mj^. 
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the lifetime of the heavy neutrino and g is the statistical 
weight of a given particle. We use gH = ,9f = 2 and 
g^ = 1, corresponding to = (/)R. The integration limits 
are 
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[EHil+4{m,/mo)y/^±PH], (8) 



where the index i = F,(j). 

For non-relativistic decays, E^ ~ n^if/2 ^ T. This 
means that the resulting fraction of low momentum 
bosons is vanishing because they are kincmatically in- 
accessible. However, for relativistic decays this is not 
the case. Here the minimum boson energy is roughly 
E^ ~ ^PH[2m'^/m'jj + ^mj^/pjj] < E^ which can be 
seen from Eq. (g). This corresponds to bosons being 
emitted backwards relative to the direction of motion 
of the parent particle. In relativistic decays it is thus 
possible to access very low energy boson states. This is 
essential for the phenomenon described below. 

Something very interesting can happen if there is ini- 
tially an overabundance of H over F. To see this, let us 
look at the structure of the decay term for the boson Eq. 
(O). Following Ref. H we make the following rough ap- 
proximation, 1 + /^ ~ Z^. The decay term then reduces 
to 



^AU calculations may be redone using g^ — 2. The changes 
are not very significant and tend to lessen the lasing effect 
described below 



Cdec [/</)] 
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(9) 



If fn > Jf this indicates an exponential growth of /<^. 
This effect has been called neutrino lasing because it is 
very similar to normal lasing, the overabundance of H 
over F corresponding to a population inversion. Exactly 
as in a laser, the process saturates when fn — fp- The 
remaining neutrinos will then decay to an almost thermal 
(j) distribution as they go non-relativistic n. 

Note that in all our actual calculations we use a zero 
initial abundance of bosons before the decay commences. 
This is not essential for the lasing phenomenon, but the 
final fraction of cold bosons is decreased if there is an 
initial thermal population (In Ref. [g| it was assumed that 
the bosons decoupled prior to the QCD phase transition 
so that their number density was depleted by a factor of 
8 relative to the other species). 

The effect will be larger for low momentum states be- 
cause of the {E^p^)~^ term. This means that if the low 
momentum states are accessible, the bosons will cluster 
in these states and produce something similar to a Bose 
condensate. For this to happen it is a necessary condi- 
tion that the decay is relativistic, because otherwise no 
low momentum states can be accessed. 

This phenomenon was first investigated in Refs. [Q-gl. 
However, Ref. Q] used a somewhat different approach, 
namely equilibrium thermodynamics, where the decays 
can lead to the formation of a true Bose condensate. In 
the above scenario, following Refs. MMi no assumption 
of equilibrium has to be made. 



III. THE 3-LEVEL LASER 

It was assumed in Refs. [p,p| that F is a hitherto un- 
known neutrino that decouples prior to the QCD phase 
transition (for example a right handed component of one 
of the known neutrinos). This is the only way of achiev- 
ing lasing in the above 2-level decay. However, it is actu- 
ally possible to obtain lasing with only the three known 
neutrinos. The effect is somewhat similar to an atomic 3- 
level laser, where particles are pumped from the ground 
state to the top level. From there they decay rapidly 
to the lasing state and from the lasing state they de- 
cay by stimulated decay to the ground state. Exactly 
the same situation can occur for neutrinos if we have a 



^Note that as mentioned we have only considered isotropic 
decays, but as mentioned in Ref. pi even if the decay is fully 
polarized and the bosons are emitted backwards, there will 
only be a rather small increase in the number of low momen- 
tum bosons. If the bosons are preferentially emitted in the 
forward direction there will be no lasing. Refs. H,h] have good 
discussions of these polarisation effects. 



three generation mass hierarchy, m^^ < m^ < m,y^ (the 
stable particles v^ and 4> are of course required to have 
masses below the cosmological limit for stable particles, 
m < 30eV, and can therefore be considered essentially 
massless in the calculations). Let us assume that we have 
the same generic decay mode as before. The possible ma- 
joron emitting decays are then n 



Vt -^ VlJi 

Vr -^ Ve 

v„ -^ Vf. 



(10) 

(11) 

(12) 



We shall assume that the second decay mode, Eq. (O), is 
suppressed so that it is negligible relative to the first, Eq. 
(|lO|). This assumption is not ruled out by the physical 
models described later. 

If the decays are non-relativistic, nothing exciting hap- 
pens and the decays are thermal. However, if all the 
decays are relativistic and the free decay of Vr proceeds 
much faster than that of v^, then we can achieve las- 
ing. The r neutrino now functions as the top level in 
a normal laser, it feeds particles to the lasing state, j/^. 
Particles in this state then decay to the ground state, v^.^ 
by stimulated decay. The condition for the lasing process 
to work is that f^ > f^^ , but because of the fast decay 
of the T neutrino this is automatically fulfilled for a large 
temperature range between rrii,^ and m^^ . 

Let us for the moment assume that the Vt decays com- 
pletely before z^^ begins to decay. It is then easy to see 
why the amount of lasing increases if Vt has decayed rel- 
ativistically. We now start with a non-thermal distri- 
bution of muon neutrinos. If Vr decays while still rela- 
tivistic, the distribution of v^ will be significantly colder 
than thermal. If we again look at the structure of Eq. 
(0), we see that for nn given, the integral increases if 
the distribution of H gets colder because of the p^ fac- 
tor in the number density integral. This means that the 
relativistic decay to muon neutrinos enhances the lasing 
from j/^ to Ve- The effect is opposite if the Vt decay is 
non-relativistic, and the lasing diminishes or disappears 
completely. 




^We choose to focus alone on the majoron decays. In princi- 
ple there might be other relevant decay modes as mentioned 
previously. 



FIG. 1. The distribution function of the light neutrino, F, 
after complete decay of the heavy neutrino, H, for different 
lifetimes of the heavy neutrino using specifically mn = 2MeV, 
mp = 0.2MeV and m^ = 0.005MeV. The fuU curve is 
for th = 0.2s, the dotted for th ~ Is and the dashed for 
TH = 20s. 

Fig. Il| shows the momentum distribution of the light 
neutrino after complete decay of the heavy neutrino. 
We see exactly that the light neutrino distribution gets 
colder if the decay is relativistic and hotter if it is non- 
relativistic. Of course, it is always a necessary condition 
for lasing to occur that i'^ decays relativistically because 
otherwise it cannot populate low momentum boson states 
at all. 



IV. NUMERICAL RESULTS 

We have investigated numerically the behavior of this 
3-level decay for different parameters by solving explicitly 
the coupled Boltzmann equations. Our numerical routine 
uses a grid in comoving momentum space and evolves for- 
ward in time using a Runge-Kutta integrator. The time 
evolution of the scale factor, R, and the photon temper- 
ature, Tj, are found by use of the energy conservation 
equation d{pR^)/dt + pd{R^)/dt = and the Friedmann 
equation H^ = SnGp/S. 
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FIG. 2. The evolution of number density for the different 
involved particles. The full curve is i^^, the dotted u^, the 
dashed Vt and the dot-dashed (j>. The parameters used are 
m^^ = 2MeV, m^^ = 0.5MeV, m^ = m^^ = 0.005MeV and 
r,, = 0.2s, r,^ = 1.0s. 

In Fig. we show an example of the number density 
evolution of the different involved particles. It is seen 
that in a large temperature interval, the Vfj, number den- 
sity is slightly higher than that of v^,. This leads to a very 
fast growth of the boson number density and is identified 
with the region where lasing takes place. The final num- 
ber density of 4> and Ve are the same because of neutrino 
number conservation (neglecting a possible small initial 
abundance of <j)) because even though the u^ abundance 
starts out higher, three 0's are produced for every two 

Ve- 

Since both v^. and cf) are almost massless there are es- 
sentially four free parameters in our equations. The first 
is the mass scale, defined for example by the r neutrino 
mass. This mass scale is unimportant for the shape of 
the spectra (this is only true if the decays take place be- 
low the weak decoupling temperature Tdcc — 2.3MeV. If 
the free decay temperature is larger than T^cc the (j) dis- 
tribution will quickly thermalize. Only well below Tjcc 
can the non-equilibrium lasing proceed, but now with an 
initial thermal (f) distribution. In our numerical work we 
choose that the free decay temperature to be lower than 
the decoupling temperature). 



FIG. 3. Final boson distribution for two different decay 
models. The full curve is for the parameters rriv^ = 2MeV, 
m^^ = 0.5MeV, m^ = m^^ = O.OOSMeV and r^^ = 1.125s, 
T^^ = 18s. The dashed curve is for m^^ = IMeV, 
m^^ = 0.25MeV, m^ = m^^ = 0.0025MeV and r^^ = 4.5s, 
T^^ — 72s. Note that although we have chosen values for m^^ 
and (jf> which is orders of magnitude larger than the experi- 
mental/cosmological limit, this does not influence the results 
as long as the values are small compared to all other relevant 
masses and temperatures. 

In Fig. H we show an example of changing the normal- 
isation {m^^) by a factor of two. This leads to an almost 
identical final boson distribution. The other three pa- 
rameters are 



3MeV 
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(13) 
(14) 
(15) 



The first two parameters describe how relativistic the free 
decays of the particles are. A thermal particle shifts from 
relativistic to non-relativistic at a temperature of roughly 
T ~ m/3. When the Universe is radiation dominated 
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Is 



IMeV/ 



Therefore if the decay is relativistic 



T < 



t(T = m/3) ~ f ^^V 



(16) 



(17) 



Thus, if Qii < 1 the decay is relativistic, whereas if a^ > 1 
it is non-relativistic. The third parameter is just the 
dimensionless v^ mass. 
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FIG. 4. The momentum spectrum of after complete 
decay for different decay parameters. Tfie full curve is 
a^^ = 0.09, a^^ = 0.01 and /3 = 0.15. The dotted is for 
a^^ — 1, civ^ = 1 and /3 = 0.25. The dashed is for a^^ = 0.5, 
a^^ = 0.5 and /3 = 0.25. The dot-dashed is for a^^ = 1, 
a^^ = 0.03 and /? = 0.25. The dot-dot-dot-dashed is for 
a^^ — 0.09, Q^^ = 1 and /3 = 0.25. The long-dashed is a 
thermal Bose distribution with the same number density. 

Fig. H shows the final boson distribution after complete 
decay for several different choices of decay parameters 
(for comparison we have also shown a thermal Bose dis- 
tribution with the same number density). It is seen that 
in order to get a large fraction of cold particles, we need 
to have ai,^ <C 1 and aj/ ^1. The reason for this was 
described in the previous section. 
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FIG. 5. Final momentum distribution of tfi. The full curve 
is for a^^ = 0.09, a^,, = 0.03 and /? = 0.25. The dotted curve 
is for a^^ = 0.09, a^^ = 0.03 and /3 = 0.15. 

Fig. H shows the effect of changing the relative mass of 
the two heavy neutrinos. We see here that to get a large 
cold fraction we also need to have rrii, /rn^^ <C 1 . This is 



also understandable if we look at the allowed momentum 
range of the produced bosons from Eq. (g) 
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T. 



(18) 



If /3 is close to one we can get very cold bosons. This 
means that the bosons will tend to cluster in low momen- 
tum states, but because of energy conservation we get a 
very hot f^( distribution. On the other hand, if /9 — > 0, 
then the lowest boson states are unaccessible and we get 
instead a colder f^ distribution. This in turn leads to 
more lasing in the lasing state and even though the frac- 
tion of cold bosons coming from the initial i>t decay gets 
smaller, this is more than compensated for by the larger 
fraction coming from the t'^ decay. This is indeed why 
the cold fraction increases when /3 decreases. 

Thus, if we want to study the region of parameter 
space with possible interest for structure formation sce- 
narios we have essentially limited ourselves to the region 
of ultrarelativistic decays with a large mass difference be- 
tween the two heavy neutrinos. It is also interesting to 
follow the evolution of the boson distributions with tem- 
perature to see how the different states populate. Fig. 

shows the temperature evolution of the boson distri- 
bution for the decay parameters a^^ — 0.09, a^ = 0.03 
and [3 = 0.25. 
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FIG. 6. The temperature evolution of the ^-distribution. 
The decay parameters are a^^ — 0.09, a,^^, — 0.03 and 
f3 = 0.25. The dotted curve is at T^ = 2.3MeV, the dashed at 
T-y = l.lMeV and the dot-dashed at T^ = 0.5MeV. The full 
curve is the final distribution after complete decay. 

The minimum possible boson energy at a given tem- 
perature is given by Eq. (Iq), where instead of ly^- we 
use v^ and for (3 we use 0. An approximate lower limit 
can be calculated by setting p^, jT^ = 3, corresponding 
roughly to the mean value from a normal thermal distri- 
bution. This lower limit is an increasing function with 



decreasing temperature and indicates roughly the mo- 
mentum where the distribution "decouples" (decoupling 
here means that the states are kinematically inaccessible) 
at a given temperature. For the model shown in Fig. |6| 
we get p7/T^ = 0.0039, 0.017, 0.081 for photon temper- 
atures of 2.3, 1.1 and 0.5 MeV. We see that the lower 
limit at a given temperature corresponds roughly to the 
momentum where the distribution decouples. Of course 
the lower limit will always be somewhat higher than the 
actual decoupling point, but it helps to understand the 
evolution of the spectra. Fig. [^ shows the corresponding 
evolution of f^. We see the before mentioned effect of 
Vt decay producing a cold Vy^ distribution and thereby 
enhancing lasing. 
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FIG. 7. The temperature evolution of the v^ distribution. 
The decay parameters are as in Fig. pi The full curve is 
the initial thermal Vp, distribution. The dotted curve is at 
T^ = 2.3MeV, the dashed at T^ = l.lMeV and the dot-dashed 
curve is at T^ = O.SMeV. 

However, even in the case where the decays are ultra- 
relativistic and the mass ratio small, the fraction of cold 
bosons is only moderate. For example in the case where 
a^^ = 0.09, Qfi/^ = 0.01 and m^^lm^^ — 0.15, the frac- 
tion of bosons with pn/T^ < 0.3 is only 44%. If the 
boson was to act as the dark matter in our Universe, a 
much larger fraction would be needed. Note, however, 
that even if the obtainable cold fraction is too small, it 
is comparable with or larger than that obtained in the 
model of Refs. 0-0] • Furthermore, it is achieved without 
invoking an unknown neutrino state. 



sufficiently large fractions of cold bosons for structure for- 
mation. 

An interesting consequence of the 3-level decay is that 
the number density of electron neutrinos is increased by 
a factor of almost three relative to the standard case. If 
the decay takes place before or during Big Bang Nucle- 
osynthesis (BBN), then it can significantly change the 
outcome because the electron neutrinos are very impor- 
tant for the neutron to proton conversion reactions. The 
higher than usual number of electron neutrinos will tend 
to keep equilibrium longer than in the standard case, 
meaning that less neutrons survive to form helium. The 
helium abundance will then be lower than in the stan- 
dard case. Furthermore, the ratio of low energy electron 
neutrinos is significantly higher than normal. This favors 
the neutron to proton reaction because of the mass differ- 
ence between the neutron and the proton and leads to the 
opposite effect of the one described by Gyuk and Turner 
|10| , where a massive r neutrino decays to electron neu- 
trinos, thereby enhancing the neutron fraction and al- 
lowing a very high baryon density without violating the 
BBN constraints. However, in their scenario, the decay- 
ing neutrino is non-relativistic and the resulting electron 
neutrinos therefore have very high energy (the opposite 
of our case) . 

Decays producing lasing effects are, however, not likely 
to take place during BBN. The reason is that neutrinos 
are kept in equilibrium by the weak reaction until roughly 
T ~ 2MeV. All neutrino distributions are then of equi- 
librium form and no lasing will take place. If the mass of 
any of the heavy neutrinos is of the same order or larger 
than this temperature, we will essentially have a non- 
relativistic decay after decoupling, but this will not pro- 
duce lasing. Thus, the heaviest neutrino in our scenario 
must be substantially lighter than 2 MeV, thus making it 
improbable that it should decay during BBN. In general 
neutrino decays affecting BBN may well take place, but 
then they will not produce lasing effects. Therefore we 
have not treated this case in the present paper. 

If we turn to more specific physical models, the most 
natural candidate for the light scalar particle, (/>, is the 
majoron. There exist several majoron models, all of them 
designed to explain the very small neutrino masses com- 
pared to the other fermions. In the simplest models, the 
decay i/2 ^ ^i + is forbidden at tree level 0, making 
the lifetime too long to be of interest. However, there 
are more complicated models where the lifetime can be 
short. In the simplest case the majoron couples to the 
neutrino via a Yukawa interaction 



V. DISCUSSION 

We have developed further the scenario originally pro- 
posed in Refs. fl-Hl by introducing a physically different 
way of obtaining neutrino lasing. However, this model 
still suffers from one of the same shortcomings that the 
original model did, namely that it is impossible to achieve 



i'-JT^ij5'^2 +h.c. 



(19) 



From this generic interaction several limits to the cou- 
pling strength, g, and therefore to the lifetime can be 
derived. For example neutrinoless double p decays give 
an upper limit on g of roughly 10~^ pT| . Other limits 
can be derived from the neutrino signal of SN1987A p2[ , 
but as mentioned in Ref. |13], it is very difficult to get a 



clear picture of the current limits. A coupling constant 
of 10"'* gives a restframe decay lifetime of [|| 
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This limit is rather weak and leaves open a large region of 
parameter space for relativistic decays in the early Uni- 
verse. Note that for our model to work the decay Vt —>■ Ve 
must be suppressed. This corresponds to the coupling gre 
being very small compared to the other coupling con- 
stants. 

Finally, we have assumed that the scattering and an- 
nihilation reactions between neutrinos and majorons are 
weak compared with the decay reaction. If, as we have 
assumed, the decay is not forbidden at tree level this is 
the case because the scattering and annihilations all are 
of order t;'*, whereas the decays are of order g^ . Were the 
annihilation and scattering reactions to be important we 
would not obtain any cold bosons because the low mo- 
mentum states would not be decoupled. Therefore, it 
is absolutely necessary that the decays are stronger than 
the scattering and annihilation reactions. Note, that even 
if the decay is of order g^ and the scattering and anni- 
hilation g^ the decay rate will be damped by a factor 
7~^ = m/ E which can be very large if the decay is ex- 
tremely relativistic. Thus, in the limit of 7 -^ cxd we have 
to account for scatterings and annihilations, but 7 would 
have to be extremely large because of the very small value 
of g. 

As far as the neutrino masses are concerned, they are 
naturally explained by see-saw models. Then only the 
lefthanded components participate in the interactions be- 
cause of the very large mass of the right handed compo- 
nents. The possible Majoron decays are then exactly 
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(21) 
(22) 



as required in our lasing model. 

In conclusion, what we have shown is that it is possi- 
ble to generate large non-equilibrium effects, even start- 
ing from identical neutrino distributions. Very few "fine- 
tunings" are necessary in order to achieve lasing. Thus, 
these effects may be present many places in early Uni- 
verse physics, both because the decay types are generic 
and because the early Universe is an inherently non- 
equilibrium environment. This was also noted in Ref. 
[|6|, but our specific model shows that the relevant range 
of models is much larger than originally thought. It is 
therefore of considerable interest to search for other phys- 
ical scenarios where the lasing process (and other similar 
non-equilibrium processes) may be effective. 
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